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3.  Report  Documentation 

3.1  Executive  summary 

The  research  on  this  project  developed  a  detailed  framework  to  describe  the  components, 
reactions,  and  internal  and  external  connections  of  cellular  control  systems.  Based  on  this 
framework,  the  principal  investigator  (PI)  built  a  stochastic  model  of  the  generalized 
biomolecular  network  and  developed  an  analytic  approach  for  the  solution.  Furthennore,  the  PI 
made  two  major  contributions  as  below: 

1.  The  PI  derived  a  novel  iterative  method  (using  both  the  forward  and  backward  iterative 
method)  to  calculate  the  solution  of  the  Network  Chemical  Master  Equation  (CME).  By 
using  a  simple  reaction  network,  the  isomerization  reaction,  the  PI  completed  the 
numerical  analysis  based  on  both  novel  analytical  method  and  direct  simulation  analysis 
for  the  test  reaction  network. 

2.  The  PI  proposed  a  novel  model  of  general  bio-molecular  networks,  where  a  number  of  different 
classes  of  the  molecules  are  involved  and  among  each  class  there  is  a  number  of  different  species 
of  molecules.  He  also  introduced  the  concept  of  the  eventual  reaction  time  (ERT)  connecting  one 
molecular  species  to  another  in  this  general  network.  Furthermore,  the  PI  derived  the 
mathematical  formulas  for  the  probability  density  function  of  the  eventual  reaction  time,  the  k-th 
moment  of  this  eventual  reaction  time,  the  probability  that  reaction  j  triggers  the  ERT  event  and 
the  joint  density  function  that  reaction  j  triggers  the  ERT  event  at  time  t. 

The  results  from  numerical  analysis  and  simulation  analysis  were  conducted  and  were  verified  to 
be  almost  exactly  same,  which  supports  the  validity  of  the  theoretical  result.  Importantly,  the 
time  to  calculate  the  solution  of  CME  and  ERT  using  our  method  is  about  10  times  faster  than 
the  time  using  the  standard  methods.  There  two  novel  methods  are  expected  to  be  superior,  in 
tenns  of  the  running  time,  for  more  complex  systems. 
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3.3  Major  Contribution 

In  all  organisms,  networks  of  coupled  biochemical  reactions  and  feedback  signals  organize 
developmental  pathways,  metabolic  pathways,  and  progression  through  the  cell  cycle.  Within 
these  regulatory  networks,  genetic  activity  is  controlled  by  molecular  signals  that  detennine 
when  and  how  often  a  given  gene  is  transcribed.  The  overall  purpose  of  this  research  project  is  to 
develop  a  new  mathematical  description  of  the  biomolecular  networks  that  control  cellular 
behavior,  formulate  both  novel  analytic  and  computational  tools  to  investigate  the  behavior  of 
these  networks  and  apply  these  approaches  to  evaluate  the  performance  of  engineered 
bioconstructs.  The  major  contribution  of  this  project  included 

>  Developed  a  general  biomolecular  network  description  of  cellular  control  systems  from 
gene  to  output  function  and  provided  a  detailed  framework  to  describe  the  components, 
reactions,  and  internal  and  external  connections  of  these  systems. 

>  Built  a  stochastic  model  of  the  generalized  biomolecular  network  and  developed  an 
analytic  approach  for  the  time-dependent  solution  and  investigate  computational  methods 
to  approximate  these  solutions. 

>  Successfully  applied  the  stochastic  model  and  analytical  methods  to  answer  specific 
questions  about  system  behavior,  including 

^  Define  the  "distance"  between  two  nodes  in  terms  of  the  time  it  takes  for  a  variation 
introduced  at  the  first  node  to  propagate  to  the  second  node.  Investigate  the  hysteresis 
effect  upon  reversal  of  the  nodes. 

•f  Compute  the  distribution  of  time  intervals  from  a  given  current  state  to  a  defined  final 
state  and  identify  the  most  likely  reaction  pathway  that  makes  this  happen. 

A  more  detailed  report  is  included  in  the  following  two  sections: 

3.3.1  A  General  Bio-Molecular  Network  Model  and  the  Solution  of  the  Network  Chemical 
Master  Equation 

The  research  of  this  research  effort  has  developed  a  modeling  framework  that  can  be  used  to 
investigate  complex  biomolecular  network  models.  Specifically,  this  research  has  reached  the 
following  goals 

>  Understood  the  molecular  basis  of  complex  stochastic  biomolecular  processes, 

>  Developed  a  stochastic  model  to  describe  these  processes  and  interactions, 

>  Derived  an  exact  solution  of  the  Chemical  Master  Equation  for  the  general  bio-molecular 
network  and  expressed  the  result  as  a  series  representation  in  terms  of  the  propensity 
functions. 

>  Simulated  biomolecular  events  during  the  operation  of  biological  constructs  under 
different  experimental  conditions. 
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At  first,  the  PI  constructed  and  defined  architecture  of  biomolecular  network  systems  as  depicted 
in  Figure  1  below. 


Figure  1.  Detailed  architecture  and  connections  between  system  components  in  our  model  of  a 

stochastic  biomolecular  network. 


Based  on  this  description,  The  PI  constructed  a  multi-dimensional  Markov  Jump  Process  on  the 
state  space 


Q  =  jv  =  (. 


-*•1,1  5  Xl,2  J  ' 


’X2,1’X2,2’' 
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in  which  i  =  1,  2,  ....  10  and  j  =  1,  2,  ....  Nr  Here  it  is  assumeed  there  are  a  total  of  M  fundamental 
chemical  reactions  in  the  network,  say  Rl,R2,---, RM  .  Each  reaction  changes  the  population  of  at  least 
one  species  of  molecules  and  a  population  change  of  any  species  comes  from  a  reaction.  Mathematically 
speaking,  reaction  Rk  is  characterized  by  a  reaction  probability  constant^,  which  is  assumed  to  be 
independent  of  time  t.  The  change  in  numbers  of  molecules  of  the  biomolecular  network  is  described  by  a 
N  =  , Nj  dimensional  vector,  called  the  state-change  vector,  vk  =  (vlk,v2k,---,vl0k)  ,  where 

vik  =  (vnk’vi  2  k’m"’ viN  k)  >  an^  Vj  j  k  is  the  change  in  the  number  of  (/, /(-molecule  produced  by 
reaction  Rk ,  for  all  i  =  1,  2,  ....  10,  j  =  1,  2,  Nt  and  k  =  1,  2,  ....  M. 


In  general,  the  propensity  functions  and  state-change  vectors  together  will  completely 
characterize  all  reaction  channels.  Mathematically  speaking,  they  also  characterize  the  stochastic 
process  X(t)  as  follows 
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P[X(t  +  dt)  =  x  +  vk\X(t)  =  x]  =  ak(x)dt  . 


(2) 


Now,  for  any  two  states  x,  y  e  Q,  we  introduce  a  new  function  of  (x,j)  as  follows: 


q 


x,y 


ak(x), 

-a0(x), 

0, 


if  there  is  a  reaction  Rk  such  that  y  =  x  +  vk , 
if  y  =  x, 
otherwise, 


(3) 


M 

where  o0(jc)  =  'Yiak(x)  for  any  x  e  Q.  Based  on  above  assumptions  for  the  reaction  rate  and 

k=l 

Equation  (1),  it  can  be  verified  that  Jf(t)  is  a  Markov  Jump  Process  and  is  defined  by  the 

10  V,. 

transition  rate  matrix  Q  =  [c/x  y  ),  with  size  of  ,9V==  ]~~[  ]~~f  1  ,  +  1 )  •  Let  P(x,  t)  be  the  probability 

;=1  7=1 

that  the  system  is  in  state  x  at  time  t,  and  P(t)  a  probability  row  vector  with  P(x, t)  as  element 
with  the  xs  ordered  lexicographically.  From  the  theory  of  stochastic  processes,  we  have 

—  P(t)  =  P(t)Q ,  with  initial  condition  P(0)  =  p0,  (4) 

dt 

where  P(0)  is  an  N  -dimensional  row  vector  detennined  by  the  initial  state  of  the  biomolecular 
network.  The  exact  solution  for  Equation  (4)  is  P(t)  =  p0eQl .  Therefore 

P(x,t)  =  P(t)ex  =  p0eQlex  ,  (5) 

where  ex  is  a  N  -dimensional  column  vector  with  one  in  the  position  of  vector  x  and  zeroes 
otherwise.  Furthermore,  we  will  have  the  following  main  result.  The  major  result  here  is  that 
the  solution  of  the  Chemical  Master  Equation  is  given  by 


P(x,t)  =  b0(x(0),x )  +  ^6„(y0),x) 


(6) 


n= 1 


n\ 


where  bn( xi0),x)  is  given  in  terms  of  the  propensity  functions  and  state  change  vectors  as: 


'/  -r ; 

[0,  otherwise. 

6„(j;<0)^)=-«o(j;<0))^-i(j;<0)^)+Zai-(j;<0))Z}»-i(j;<0)  +vk’x)’  for  n  =  1,2,  ••• 


(7) 


k—\ 


where  J(0)  is  any  state  in  the  state  space  Q. 
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3.3.2.  Eventual  Reaction  Time  in  General  Bio-Molecular  Networks 

The  PI  proposed  a  novel  model  of  general  bio-molecular  networks,  where  a  number  of  different 
classes  of  the  molecules  are  involved  and  among  each  class  there  is  a  number  of  different  species 
of  molecules.  The  major  contribution  in  this  part  included: 

>  Introduced  the  concept  of  the  eventual  reaction  time  (ERT)  connecting  one  molecular 
species  to  another  in  this  general  network; 

>  Derived  the  mathematical  fonnulas  for  the  probability  density  function  of  the  eventual 
reaction  time,  the  k-th  moment  of  this  eventual  reaction  time,  the  probability  that  reaction 
j  triggers  the  ERT  event  and  the  joint  density  function  that  reaction  j  triggers  the  ERT 
event  at  time  t. 

At  first,  the  PI  considered  this  ERT  problem  from  the  viewpoint  of  the  cell  network.  Since  there 
are  many  other  types  and  classes  of  the  molecules  at  time  t0 ,  we  may  assume  that  at  t()  there  are 

also  a  specific  number  xPoq0(to)  °f  (po,qo)-molecules,  in  which 
Po  =1,2,  •••,10;  q0  =  1,  2,  ■■■,  N%  where  (p0  ,q0)^  (/0 ,  j0 )  .  The  set  of  all  of  these  values  taken 
together  describe  the  state  of  the  system  at  time  t0.  We  denote  by  x0  this  network  state  at  time 
t0 ,  which  is  a  specific  state  of  the  state  space  Q  (Equation  (3)),  and  is  called  the  initial  state  of 
the  network.  Since  there  are  many  other  species  from  different  classes  of  the  molecules  at  any 
specific  time,  at  the  time  when  the  system  first  attains  a  specific  number  x,.  ;  of  (iiji)- 

molecules  in  the  network,  it  is  possible  to  have  a  number  x  of  (p/.^yj-molecules,  in  which 
Pi  =1,2,  •••,10;  qx  =1,  2,  •••,  Nqi  and  (/?,  ^  (/,,/, )  •  Any  of  these  states  could  describe  a 

legitimate  destination  state.  Since,  we  are  only  interested  in  the  first  time  the  system  contains 
exactly  x.  .  (//.//(-molecules,  we  may  form  a  sub-set  of  the  state  space,  ClD ,  consisting  of  all 

those  states.  For  example,  if  (/, ,  /j )  =  (1,2) ,  then 


Q  —  (xj  j ,  x.  ,  •  •  • ,  ,  x2l ,  x22 ,  •  •  • ,  x2 ,  •  •  • ,  x10  j ,  . 


10,2’ 


JeQ}- 


This  set  is  called  the  set  of  destination  states.  If  we  denote  by  ClE  the  set  of  all  other  state,  then 
Q  =  Q£+Qd  .  Accordingly,  we  use  cf,  af,  RE ,  ME  to  denote  the  reaction  probability 
constant  for  reaction  k,  the  corresponding  propensity  function,  the  reactions  and  the  total  reaction 
numbers  in  ClE .  Based  on  above  two  steps  description,  if  we  denote  by  r(  ^  the  time  interval 

when  the  system  first  attains  the  specific  number  x(.  ■  of  (//.//(-molecules,  then  what  is  the 
probability  density  function,  Pix  D.(r),  of  this  random  variable  r,x  _^D, ,  and  what  is  the  k-th 


moment  of  z, 


Uo 


? 

>D)  ' 


If  we  define  the  set  of  states  D  as  absorbing  states,  which  means  that  the  system  dynamics  will 
stop  as  soon  as  the  system  reaches  any  member  of  the  set  D,  then  the  system  stopping  time, 
starting  from  the  state  x0 ,  will  be  _>D) .  We  will  briefly  introduce  our  methods  and  the  results. 
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First,  we  introduce  a  new  matrix  S,  which  is  obtained  from  the  matrix  Q  by  deleting  all  rows  and 
columns  corresponding  to  states  in  D.  Secondly,  we  need  to  detennine  the  number  of  states  in  D, 
denoted  Dn,  and  introduce  an  ( 9VDn ) -d  i  m en  s i  o n a  1  row  vector  6 \  which  is  1  in  the  corresponding 
position  of  x0  and  zero  otherwise.  Based  on  the  concept  of  phase  type  distribution  introduced  by 
Neuts  (1982),  we  can  derive  the  probability  distribution  of  the  eventual  reaction  time  (ERT)  by 
p(r(  _^D)  <  t)=  1  -9es'e ,  t  >  0,  where  e  is  a  (,9fD»)-dimensional  column  vector  with  one  in  all 
elements. 


Furthermore,  by  several  mathematical  arguments,  we  may  now  derive  and  summarize  our  main 
results  as  follows: 

>  The  probability  density  function  of  the  eventual  reaction  time  (ERT)  is 

oo  sit-l 

/>(,„./» (0  =  -Oes,Se  =  '  1  -  °> 

where  e  is  a  (N-Dn)-dimensional  column  vector  with  one  in  all  elements,  and  b^(x0 )  is 
given  in  terms  of  the  propensity  functions  and  state  change  vectors  as  follows: 

me  ud 

bi(xo)  =  -ao(xo)+  Tjak(xo)=  -«oB(^)  =  -Z  ak(xo) 

k= 1  k=\ 

me 

K ix0 )=~ao(xo) bL (*r)+  Z ak  ( xo ) K- 1  (*0  +  vk ) »  for  n  =  2, •  •  •  . 

k= 1 

>  The  A-th  moment  of  this  eventual  reaction  time  (ERT)  is: 

E[t(Xo^D)T  =  (-\)kk\(dS  ke),  k=  1,2,3,  .... 

>  If  we  further  ask  “which  reaction  triggers  this  ERT  event?”  we  have  the  following  result. 
Let  R i  be  the  set  of  reactions  such  that  after  each  reaction  of  if;  a  (//,//)-moleculc  is 
produced  as  a  product.  Therefore, 

r(  reaction  Rj  triggers  this  ERT  event)  = - - - ,  R/  e  Rt . 

Z/teif(1>  Ck 


> 


If  we  denote  by  „,(/,/)  the  joint  density  function  that  the  eventual  reaction  time  is  t 

and  results  as  a  consequence  of  reaction  j,  i.e., 

p(f  <  r(v>/)|  <t  +  dt  and  reaction  R;  triggers  this  event)  =  P{XoD)  (t,  j)dt 

r  \ 


then 


(tj)  = 


I 


2X(*.) 


n=\ 


in- 1)! 


for  t  >  0. 
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3.3.3.  Case  Study. 

By  using  a  simple  reaction  network,  the  isomerization  reaction,  The  PI  reached  the  following 
conclusion. 

•  Completed  the  numerical  analysis  based  on  both  our  novel  analytical  method  and  direct 
simulation  analysis  for  the  test  reaction  network. 

•  The  results  from  numerical  analysis  and  simulation  analysis  are  almost  exactly  same, 
which  supports  the  validity  of  our  theoretical  result. 

•  The  time  to  calculate  the  solution  of  CME  using  our  method  is  about  10  times  faster  than 
the  time  using  the  standard  methods. 

•  Our  method  is  expected  to  be  superior,  in  terms  of  the  running  time,  for  more  complex 
systems.  This  is  currently  being  tested. 
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3.4  Matlab  Code  and  the  corresponding  diagrams 


A.  The  Matlab  code  for  calculating  the  Network  Chemical  Master  Equation  (CME)  by  using 
the  new  proposed  iterative  method  is  as  follows: 


function  PxT  =  GetCMEPxT() 

%GETPXTS  Summary  of  this  function  goes  here 
%  Detailed  explanation  goes  here 
C  =  .2; 

X0  =  (0:1:10); 

XI  =(10:-1:0); 

A1  =  X0*C; 

A2  =  X1*C; 

AO  =  A1+A2; 

N=  1000; 
t  =  (0: .  1 : 1 0); 

curSta  =  0;  %Initial  State  PcurSta(x,t) 

PxT  =  zeros(length(X0),  length(t)); 
for  tV  =  1 :  length(t) 
for  x  =  1  :  length(XO) 

%Get  SumD 

SumD  =  GetDn(N,  x,  Al,  A2,  A0,  t(tV)); 

PxT(x,tV)  =  SumD(N,  curSta+1);  %Our  Index  start  at  1  not  0  so  we  add  1 
end 
end 

plot(t,  PxT); 
xlabel(’Time  t'); 

ylabel(sprintf('P%d(x,t)',  curSta)); 
title(sprintf('P%d(x,t)',  curSta)); 

legend(sprmtf('P%d(X0,t)', curSta),  sprintf('P%d(X  1  ,t)', curSta), 
sprintf('P%d(X2,t)', curSta),  sprintf('P%d(X3,t)', curSta),  sprintf('P%d(X4,t)', curSta), 
sprintf('P%d(X5,t)', curSta),  sprintf('P%d(X6,t)', curSta),  sprintf('P%d(X7,t)', curSta), 
sprintf('P%d(X8,t)', curSta),  sprintf('P%d(X9,t)', curSta),  sprintf('P%d(X  1 0,t)', curSta)); 

%filNam  =  sprintf('Result/PxT.xls'); 

%celNam  =  sprintf('PxT’); 

%xlswrite(  filNam,  PxT,  celNam,  Al') 

function  SumD  =  GetDn(N,  y,  Al,  A2,  A0,  t) 

MaxStates  =  length(AO); 

D  =  zeros(N,  MaxStates); 
for  n  =  1 :  N 

for  x  =  1 :  MaxStates 
if  n  ==  1 
if  x  ==  y 

_ D(n,x)=  1; _ 


11 


Final  Performance  Report 


else 

D(n,x)  =  0; 
end 
else 

D(n,x)=  -A0(x)*D(n-l,x); 
if  x  >  1 

D(n,x)  =  D(n,x)  +  Al(x)*D(n-l,  x-1); 
end 

if  x  <  MaxStates 

D(n,x)  =  D(n,x)  +  A2(x)*D(n-l,  x+1); 
end 

D(n,x)  =  D(n,x)  *  t  /  n; 
end 
end 
end 

SumD  =  zeros(N,  MaxStates); 
for  n  =  1 :  N 

for  x  =  1 :  MaxStates 
if  n  ==  1 

SumD(n,x)  =  0; 
else 

SumD(n,x)  =  SumD(n-l,x)  +  D(n,x); 
end 
end 
end 

for  n  =  1 :  N 

for  x  =  1 :  MaxStates 
if  y  ==  x 

SumD(n,x)  =  SumD(n,x)  +  1 ; 
end 
end 
end 


(l'x)Od 
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The  corresponding  diagrams  are: 


PO(x.t) 
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B.  The  Matlab  code  for  calculating  the  network  Eventual  Reaction  Time  (ERT)  by  using  the 
new  proposed  iterative  method  is  as  follows: 


function  PxT  =  GetERTPxT() 

%GETPXTS  Summary  of  this  function  goes  here 
%  Detailed  explanation  goes  here 
Cl  =  .2; 

C2  =  0.05; 

X0  =  (0:1:10); 

XI  =(10:-1:0); 

A1  =X0*C1; 

A2  =  X1*C2; 

AO  =  A1+A2; 

N  =  450; 
t  =  (0:1:100); 

curSta  =  0;  %Initial  State  PcurSta(x,t) 

PxT  =  zeros(length(X0),  length(t)); 
for  tV  =  1 :  length(t) 
for  x  =  1  :  length(XO) 

%Get  SumD 

SumD  =  GetCn(N,  x,  Al,  A2,  A0,  t(tV)); 

%PxT(x,tV)  =  SumD(N,  curSta+1);  %Our  Index  start  at  1  not  0  so  we  add  1 
end 
end 

%  plot(t,  PxT); 
plot(SumD); 

%  xlabel(’Time  t'); 

%  ylabel(sprintf('P%d(x,t)',  curSta)); 

%  title(sprintf('P%d(x,t)',  curSta)); 

%  legend(sprintf('P%d(X0,t)', curSta),  sprintf('P%d(X  1  ,t)', curSta), 
sprintf('P%d(X2,t)', curSta),  sprintf('P%d(X3,t)', curSta),  sprintf('P%d(X4,t)', curSta), 
sprintf('P%d(X5,t)', curSta),  sprintf('P%d(X6,t)', curSta),  sprintf('P%d(X7,t)', curSta), 
sprintf('P%d(X8,t)', curSta),  sprintf('P%d(X9,t)', curSta),  sprintf('P%d(X10,t)', curSta)); 

%filNam  =  sprintf('Result/PxT.xls'); 

%celNam  =  sprintf('PxT'); 

%xlswrite(  filNam,  SumD,  celNam,  ’AT) 

function  SumC  =  GetCn(N,  y,  Al,  A2,  A0,  t) 

MaxStates  =  length(AO); 

C  =  zeros(N,  MaxStates); 
for  n  =  1 :  N 

for  x  =  1 :  MaxStates 
if  n  ==  1 
if  x  ==  9 

_ C(n,x)  =  t  *  (-A2(x)); _ 
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else 

C(n,x)  =  0; 
end 
else 

C(n,x)=  -A0(x)*C(n-l,x); 
if  x  >  1 

C(n,x)  =  C(n,x)  +  Al(x)*C(n-l,  x-1); 
end 

if  x  <  MaxStates 
if  x  ~=  9 

C(n,x)  =  C(n,x)  +  A2(x)*C(n-l,  x+1); 
end 
end 

C(n,x)  =  C(n,x)  *  t  /  n; 
end 
end 
end 

SumC  =  zeros(N,  MaxStates); 
for  n  =  1 :  N 

for  x  =  1 :  MaxStates 
if  n  ==  1 

SumC(n,x)  =  0; 
else 

SumC(n,x)  =  SumC(n- 1  ,x)  +  C(n,x); 
end 
end 
end 

% 

%for  n  =  1 :  N 
%  for  x  =  1 :  MaxStates 
%  if  y  ==  x 

%  SumC(n,x)  =  SumC(n,x)  +  1; 

%  end 

%  end 
%end 
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The  corresponding  diagrams  are: 
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4.  Conclusion 

The  research  of  this  project  developed  a  new  mathematical  description  of  the  biomolecular 
networks  that  control  cellular  behavior,  formulate  both  novel  analytic  and  computational  tools  to 
investigate  the  behavior  of  these  networks  and  apply  these  approaches  to  evaluate  the 
performance  of  engineered  bio-constructs.  The  PI  has  investigated  and  developed  a  novel  general 
stochastic  biomolecular  network  description  of  cellular  control  systems  from  gene  to  output 
function.  He  provided  a  detailed  framework  to  describe  the  components,  reactions,  and  internal 
and  external  connections  of  these  systems,  built  a  probability  model  of  the  generalized  stochastic 
biomolecular  network,  developed  an  analytic  approach  for  the  time-dependent  solution  and 
investigated  computational  methods  to  approximate  these  solutions.  Specifically,  the  PI 
successfully  applied  the  stochastic  model  and  analytical  methods  to  have  answered  specific 
questions  about  system  behavior,  which  include  define  the  "distance"  between  two  nodes  in 
terms  of  the  time  it  takes  for  a  variation  introduced  at  the  first  node  to  propagate  to  the  second 
node,  investigated  the  hysteresis  effect  upon  reversal  of  the  nodes,  and  computed  the  distribution 
of  time  intervals  from  a  given  current  state  to  a  defined  final  state  and  identify  the  most  likely 
reaction  pathway  that  makes  this  happen.  By  successful  dissemination  of  these  technical 
advances,  further  technical  progresses  can  be  expected  by  collaborations  and  technical 
interactions  with  colleagues. 

The  project  partially  supported  two  Master  students  in  the  Department  of  Computer  Science  at 
Texas  Southern  University,  one  of  the  major  Historically  Black  Colleges  and  Universities 
(HBCU)  in  the  nation  with  high  concentration  of  minority  students  from  under-represented 
groups.  The  project  directly  impacted  the  capabilities,  interests,  and  careers  of  the  graduate 
students  who  participate  in  the  project  as  well  as  any  other  students  of  interest. 
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